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Abstract We propose a three-stage model with Blandford- 
Znajek (BZ) and hyperaccretion process to interpret the re- 
cent observations of early afterglows of Gamma-Ray Bursts 
(GRBs). In the first stage, the prompt GRB is powered by a 
rotating black hole (BH) invoking the BZ process. The sec- 
ond stage is a quiet stage, in which the BZ process is shut 
off, and the accretion onto the BH is depressed by the torque 
exerted by the magnetic coupling (MC) process. Part of the 
rotational energy transported by the MC process from the 
BH is stored in the disk as magnetic energy. In the third 
stage, the MC process is shut off when the magnetic energy 
in the disk accumulates and triggers the magnetic instabil- 
ity. At this moment, the hyperaccretion process may onset, 
and the jet launched in this restarted central engine gener- 
ates the observed X-ray flares. This model can account for 
energies and timescales of GRBs with X-ray flares observed 
in early afterglows. 
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1 INTRODUCTION 



Recently, the flares in the early X-ray afterglows have been discovered in 
both long and short bursts with the X-ray telescope (XRT) on board Swift. 
These flares appear at about tens to hundreds seconds after the trigger of 
the GRBs and last about several hundred seconds (for review of the Swift 
results and, in paricular, the X-ray flares, see Zhang 2007). Their rapid 
variability has been interpreted as the inner engine being active much 
longer than the duration of the GRB itself (Zhang et al. 2006). Therefore, 
it is needed that after the cease of the prompt gamma-ray emission, the 
central engine can be restarted (Fan & Wei 2005; Zhang et al. 2006). 
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For the case of long GRBs, King et al. (2005) suggested that the X- 
ray flares could be produced from the fragmentation of the collapsing 
stellar core in a modified hypernova scenario. The fragment subsequently 
merges with the main compact object formed in the collapse, releasing 
extra energy. In this two-stage collapse model, the time delay between the 
burst and the flare reflects the gravitational radiation time scale for the 
orbiting fragment to be dragged in. 

For the case of short GRBs, MacFadyen, Ramirez-Ruiz & Zhang (2005) 
suggested that the flares could be the result of the interaction between the 
GRB outflow and a non-compact stellar companion, in this model, short 
GRBs result from the collapse of a rapidly rotating neutron star in a close 
binary system. While Dai et al. (2006) suggested that the flares in short 
GRBs can be produced by differentially rotating, millisecond massive pul- 
sars after the mergers of binary neutron star. The differential rotation leads 
to windup of interior poloidal magnetic fields and the resulting toroidal 
fields are strong to float up and break through the stellar surface. Magnetic 
reconnection-driven explosive events then occur, leading to multiple X-ray 
flares minutes after the original gamma-ray burst. Gao & Fan (2006) also 
suggested a short-lived supcrmassive magnetar model to account for the 
X-ray flares following short GRBs. In their model, the X-ray flares are 
powered by the dipole radiation of the magnetar. 

Motivated by the fact that the flares are observed in both long and 
short duration GRBs, Perna et al. (2005) suggested that the flares could 
be produced due to subsequent accretion of blobs of material in the hyper- 
accreting accretion disk, which initially circularize at various radii and sub- 
sequently evolve viscously. While Proga & Zhang (2006) conjectured that 
the energy release can be repeatedly stopped and then restarted by the 
magnetic flux accumulated around the accretor. The restarting of the ac- 
cretion corresponds to the X-ray flares. 

Although the main features of the X-ray flares can be interpreted suc- 
cessfully in the above scenarios, a quantitative model with a clear physical 
process for restarting the central engine remains lacking. 

Very recently, Lei et al. (2005a, hereafter L05; 2005b) proposed a 
scenario for GRBs in Type Ib/c SNe, invoking the coexistence of the 
Blandford-Znajek (BZ; Blandford & Znajek 1977) and Magnetic Coupling 
(MC; Blandford 1999; van Putten 1999; Li 2000, 2002; Wang et al. 2002) 
processes. The BZ process can provide "clean" energy for the prompt GRB 
by extracting rotating energy from a BH (Lee et al. 2000). In L05, the ac- 
cretion onto the BH is depressed by the torque exerted by the BH through 
MC process. It is found that the BZ process will be shut off before the MC 
process because of the close of the half-opening angle on the BH horizon. 
Recently, van Putten & Levinson (2003, hereafter PL03) suggested that 
the magnetic energy accumulated by the MC process can trigger the mag- 
netic instability, which results in the cease of the MC process (Eikenberry 
& van Putten 2003). 

Based on L05 and PL03, we suggest that the pause and the restart 
of the central engine arise from the MC process, and propose a three- 
stage model to interpret the GRBs and their flares. In this model, the 
prompt GRB and the X-ray flare are powered by the BZ process and the 
subsequent hyperaccretion, respectively. 
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Fig. 1 A schematic drawing of the magnetic field configuration of 
a three-stage model 



2 THE THREE-STAGE MODEL 

In this section, we present the three-stage model for the inner engine for 
GRBs. The magnetic field configuration is shown in Figure 1, which is 
adapted from van Putten (2001). 

In Figure 1 the angle 9h is the half-opening angle of the open magnetic 
flux tube, indicating the angular boundary between open and closed field 
lines on the horizon. The angle Oh can be determined by (Wang et al. 
2003) 

cosOh = G(a^tn)dC (1) 

where, a» = J/M 2 is the BH spin defined in terms of the BH mass M 
and the angular momentum J, the parameter n is the power-law index for 
the variation of B£, i.e.,B£ oc £~™,and £ = r/r ms is the radial coordinate 
on the disk, which is defined in terms of the radius r ms = Mxms °f the 
marginally stable orbit (Novikov & Thorne 1973). Xms depends only on 
(I* as (Novikov & Thorne 1973) 

Xms ~ 6 xL + 8a*Xm S - 3a 2 = 0. (2) 
The function G (a*; £, n) is given by 

e- n X 2 ms V 1 + «*X™4" 2 + 2a2 X ™ 6 4-3 
G(a*;£,n) = — ^=. 

2^(1 + alx m i + 2alxVn%) (1 - 2 Xm 2 s £-i + alx^ 2 ) 

(3) 

The parameter £s = rs/r ms in the equation |T]),is the critical radius 
of screw instability in the MC process, which is determined by the criteria 
for the screw instability given by Wang et al. (2004), i.e., 

^ < ^ (4) 
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Fig. 2 The critical lines for screw instability in parameter space, 
by which the space is divided into regions I and II. Each filled 
circle with arrowhead in the regions represents one evolutionary 
state of the BH. 



where L is the poloidal length of the closed field line, and vjd is the 
cylindrical radius on the disk, B£ and Bp are the poloidal and toroidal 
componets of the magnetic field on the disk, respectively. 2ttzud/L and 
Bp I Bp are expressed as (Wang et al. 2004) 



BllB\ 



2a* (1 - (3) 



^-"(l + q)[2csc 2 6»-(l- 9 )] 1 



r 2 + 2alx™U~ 3 



1 + alxms + 2alx 



6 
rn.s 



(5) 



2-KU]£) 



alt, 2 Xm 4 s 



(6) 

By using equation (TTJ) we have the critical lines for the half-opening 
angle Oh = in a* — n parameter space as shown in Figure 2. 

An interesting feature shown in Figure 2 is that the angle Oh can 
evolve to zero (the shaded region) with the decreasing a* for 2.616 < n < 
5.493. This evolution characteristic implies that the BZ process will be shut 
off when the BH spin decreases to the critical value a^ RB corresponding 
to Oh = 0. The lifetime of the half-opening angle Oh is defined as the 
evolution time of BH from the initial spin a* (0) to af RB . 



2.1 Stage 1: the Prompt Emission of GRB 

In the first stage, the BZ process works together with the MC process and 
the evolution state of the BH is located in region I of Figure 2. 

Since angular momentum is transferred from the rapidly rotating BH 
to the disk, a very strong torque will be exerted on the inner disk in the 
MC process. The accretion onto the BH will be depressed by the MC 
process (van Putten & Ostriker 2001, PL03; Li 2002). Therefore, in this 
stage, the accretion rate is very low, and the prompt emission of GRB is 
mainly powered by the BZ process. 
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The powers and torques for the BZ and MC processes are expressed 
as (Wang et al. 2002; Wang et al. 2003) 



)h k (1 - k) sin 3 
2-(l-q) sin 2 



Pbz = Pbz/Po = 2al / ^ M ' (?) 



T BZ = T BZ /T a = 4a*(l + q) I £ ^ ' - - (8) 



' H (l-fc)sin 3 * 
2 - (1 - (?) sin 2 9 1 



-"."7 • / 2 _ (1 _ a)sin ^< 

T T "r in • ^ (I-/?) sin 3 gdg 

T A/C = Tmc/To = 4a* (1 + g) / - — . 2 , (10) 

V H 2 - (1 - g) sin 



where we have g = y/l - of, P = PfjAf 2 and T ee B 2 h M 3 . B h is the 
magnetic field on the BH horizon. The parameters k and /3 are the ratios 
of the angular velocities of the open and closed magnetic field lines to 
that of the BH, respectively. Usually, k — 0.5 is taken for the optimal BZ 
power. 

The energy E B z extracted in the BZ process can be calculated during 
the evolution of GRB. The true energy for GRBs, P 7 , is given in L05, 



£L = e~E B z = 2.69xlO D3 erasx — y- —L- / — = da* 

7 7 V M Q ) \ 0.15; 1.(0) -T mas + 2a*P mas 

(11) 

where P mag ee P bz + P Ma , f mog ee f BZ + f MC , and M ee M/M(0) is 
the ratio of the BH mass to its initial value. Following van Putten et al. 
(2004), we take the efficiency e 7 = 0.15 in calculations. 

The duration of GRB, tcRB, is defined as the lifetime of the half- 
opening angle 0h, which is exactly equal to the time for the BH spin 
evolving from a*(0) to a* B , i.e., 



t G RB = 2.7 x 10 3 s x — ^- / — = = da*, 

\ B H ) \M{0)) J am -T mag + 2a*P mag 

(12) 

The energy and duration of the prompt emission can be obtained 
by taking M(0), a*(0), n and P# into equations (fTTj) and (TT2"|) . In the 
following calculations the initial BH mass and spin are assumed to be 
M(0) = 7M Q and a*(0) = 0.9. 

By adopting proper values of parameters ro and P#, our model can ex- 
plain both long and short bursts, such as GRB 050219A and GRB 050709 
(see Table 1). 

From Table 1 we find that the shorter duration and lesser energy of 
short burst could be fitted with smaller value of n and greater value of 
Bh- Some long bursts have been fitted in L05 by reasonable values of n 
and Bh- 
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Table 1 GRB 050219A and GRB 050709. 



GRBs 


Type 


T 90 (s) 


£ 7 (10 50 er ff ) 


n 


Bh(10 15 G) 


50219A 


Long 


23.6" 


> 4.8° 


>3.545 


>0.7 


050709 


Short — hard 


0.073 


0.021 6 


3.330 


2 



a Tagliaferri et al. 2005 
b Fox et al. 2005 



2.2 Stage 2: Accretion Depressed by the MC Process 

In this stage, the MC process is still active, and it depresses the accretion 
onto the BH. Following PL03, the magnetic field energy of disk eb will 
build up in response to the power received from the BH, i.e., 

eu(t)= / vPAicdt, (13) 
Jo 

where, r\ is the ratio of the increase rate of the magnetic field energy to 
the total MC power, and 77 <gc 1 since most of the rotational energy of a 
Kerr BH is emitted in gravitational radiation and MeV neutrino emissions 
(van Putten 2001). 

According to PL03 the magnetic instability of the disk will occur if 

e B /e K > 1/15, (14) 

where Ek is the kinetic energy of disk. Once this instability is triggered, the 
magnetic field lines connecting the BH horizon and the disk will disconnect 
quickly, i.e., the MC process will be ceased. 

As the accretion in the inner part of the disk is blocked by the MC 
process the inner region will pile up at the accretion rate corresponding 
to the radius rs, and a torus might form due to the matter accumulation. 

In the case without the depression of the MC process the accretion rate 
can be estimated by considering the balance between the pressure of the 
magnetic field and the ram pressure of the innermost part of the accretion 
flow (Moderskin et al. 1997), i.e., 

B 2 h /(8tt) = P ram ~ pc 2 ~ M D /{A-KT 2 H ), (15) 

As seen in Figure 1, the MC torque mainly affects the accretion flow 
within the radius £s. As a simple model, we assume the accretion rate of 
the disk beyond £s can also be estimated by equation ([15"]) . The time of 
forming torus in the suspended accretion state can be estimated as 

Ttorus pa = aB^mg 1 - -j2.7 x 10 3 s, (16) 

M D (l + <?) 

where Mr is the mass of torus, and run = M/Mq, a = Mt/M, 
B 15 = B H /10 15 G. For M T = O.1M , B H = 10 15 G and a* = 0.5, we 
have Ttorus ~ 3.2s. Thus a torus will form after (for a short bursts) or 
during (for a long bursts) the prompt emission of GRB. The kinetic en- 
ergy of the torus Ek can be estimated as (PL03) 



_ M T M 
~ 2r T ' 



(17) 
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where, the radius of the torus tt can be estimated by ~ r ms . 

Once the magnetic field energy of disk eb satisfies equation (| 14[) . the 
instability occurs, and the magnetic field lines connecting the BH horizon 
and the disk will be disconnected quickly. At this time, the MC process 
ceases, and the BH spin is denoted as a° cc . 

2.3 Hyperaccretion as the Restarting Engine for X-Ray Flares 

The torque to depress the accretion flow will disappear immediately once 
the MC process stops. At this stage the viscous torque is dominative, and 
the hyperaccretion will occurs due to an abrupt lose of angular moment 
of disk matter. 

Since the BZ process has been ceased at the end of the first stage due 
to the shutting off of the half-opening angle Oh, the only operating energy 
process at this time (the third stage) is hyperaccretion. 

We assume that the X-ray flares are powered by the hyperaccretion, of 
which the starting time is denoted as t acc (in the frame of central engine). 
In addition, we assume that the starting time t acc is determined by the 
duration for the BH spin evolving from a*(0) to a" cc . 

The observational beginning time of the flare is tfi ar e -We relate t aC c 
to the radius for the new ejecta (produced by the hyperaccretion process) 
to catch up with the foregoing merged shell of GRB prompt stage by 



Here the Lorentz factors of the ejected shells and their diversity are 
all assumed in the same order 7(Sari & Piran 1995). Considering that the 
observational time obeys tfi are ss R/j 2 c, we can regard t acc as the time 
of observing the X-ray flare, i.e., tf t 

are ^ ^acc- 

In our model t acc is sensitive to the values of n, Bh and rj, where n and 
Bh are obtained by fitting the energy and duration of the GRB prompt 
emssion (See the discussion in stage 1). Therefore, r\ is only determined by 
tacc or t f i are for a given GRB. For GRB 050219A and GRB 050709, the 
X-ray flares are observed at i/; are ~ 100s after the trigger of the bursts. 
The values of the parameter are taken as r] = 0.05 and 0.01 for GRB 
050219A and GRB 050709 in calculations, respectively. 

For other bursts, the values of rj are in the same order. This can be 
proved by the following estimation. For Mt = O.IMq, the critical mag- 
netic energy is Eb = « 10 51 erg. The MC power can be estimated 
as P M c ~ B 2 M 2 Rj 6.59 x I0 50 erg/s '. lit nares ~ 100s, we have rj ~ 0.02. 

It has been argued that the hyperaccreting torus can power the X-ray 
flares via neutrino annihilation (e.g. Ruffert & Janka 1998). The jet lumi- 
nosity driven by this mechanism depends on the mass accretion rate very 
sensitively, since the neutrino emission is related very closely to the den- 
sity and the temperature of the torus. For accretion rates (Md) between 
0.01 and 0.1M e s -1 , the vv annihilation luminosity can be well fitted by 
(Popham et al. 1999; Fryer et al. 1999; Janiuk et al. 2004) 



R f ~ 7 2 ci, 



(18) 




(19) 



In this stage, the BH spin must be aj cc , which is about 0.5. 
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The total energy of the flare is also typically smaller than that of 
the prompt emission, although in some cases both could be comparable 
(e.g. for GRB 050502B). Based on equation (fl"9|) . a torus with mass of 
O.IMq can provide 10 50 erg via neutrino annihilation if its accretion rate 
is about O.IMq/s. This energy is enough to power the X-ray flares in 
GRB afterglow. 

The duration of the flare can be estimated as following. As suggested 
by Zou, Dai & Xu (2006), we assume that the emission of X-ray flares come 
from the reverse shock of internal shock. For a relativistic reverse shock, 
the duration of the shock in the observer's frame is Atfi are ~ A/c(Sari & 
Piran 1995), while the depth of shell for the spreading case is A ~ R/l 2 - 
Thus, the duration of one flare is 

Atfiare ^ Rf/j 2 C ~ t flare (20) 

It is interesting to note that the three time scales are equal approxi- 
mately, Atfi are ~ t fiare ^ tacc: and they are all independent of the Lorentz 
factors. This result is also in agreement with the observations. 

Usually, the afterglow light curves are characterized by multiple flares. 
This feature could be understood based on restarting the hyperaccretion 
several times, provided that the disk may be fragmentized by the gravita- 
tional instability as proposed by Perna et al. (2005), or the accretion can 
be repeatedly stopped and then restarted by the magnetic flux accumu- 
lated as shown by Proga & Zhang (2006). 

Moreover, the durations of these flares seem to be positively correlated 
with the epochs when the flares happen, i.e., the later the epoch, the longer 
the duration. There is particularly evident in GRB 050502B and the short- 
hard GRB 050724. In both cases, there is an early flare (several 100s for 
GRB 050502B and several 10s for GRB 050724) whose duration is of the 
order of the peak time itself, and there is also a very late flare at tens of 
thousands seconds with a duration of the same order. The above properties 
of the flares could be interpreted by equation (f2T)]) . which predicts that 
the arrival time of each new flare should directly correlate with its total 
duration. 

3 DISCUSSION 

In this paper, we propose a three-stage model for the inner engine for 
GRBs. The prompt emission of GRB is explained by the BZ process, and 
the X-ray flares are explained by the latter hyperaccretion process. One 
of the features of our model is the MC effects on the accretion process, 
both the depress accretion and the start of the hyperaccretion process are 
related intimately to the MC mechanism. It turns out that some features 
of the X-ray flares in the early afterglow after prompt gamma ray emission 
can be interpreted. 

Comparing with the models of Perna et al. (2005) and Proga & Zhang 
(2006), our model can describe the cease and the restart of the central 
engine in a quantitative way. While the former have their advantages in 
explaining the multiple flares. Our model can also interpret the multiple 
flares if the MC process can be stopped and restarted repeatedly. Although 
the restart of the MC process is very speculative, it is possible since the 
magnetic filed can be recovered by the dynamo process (Hawley et al. 1995) 
or the accumulation of the magnetic flux (Spruit, Stehle and Papaloizou 
1995). 
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In our model, the X-ray flares are powered by the hyperaccretion pro- 
cess. However, for short GRBs, Fan et al. (2005) suggested that this pro- 
cess is not efficient. Instead, the BZ process is a very powerful mechanism. 
Similar results are found in Xie et al. (2007). Like the MC process, the 
restart of the BZ process is also possible. 

From the above discussion we find that the rebuilding of the magnetic 
field is essential to our model, and we intend to address this issue in our 
future work. 
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